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<C ! Abstract 



We present an exactly solvable model for photon emission, which allows us 
to examine the evolution of the photon wave function in space and time. We 



> 
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■ apply this model to coherent phenomena in three level sytems with special 

emphasis on the effects of the photon detection process. 

C^. I. INTRODUCTION 
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Describing the open system dynamics of a quantum object, the states of external fields 



such as the light field are usually represented by some form of statistical approximation: for 
instance, the master equation approach assumes an unchanging thermal state of the external 
fields, disregarding all effects which the system might have on them ("bath approximation"). 

Approaches such as the Monte Carlo wave function method or the subensemble density 
matrix consider measurements of the photons at the instant of their emission, thereby 
arriving at an easy to handle description of the open system dynamics while avoiding a more 
detailed description of the light field itself. 

However, in typical experimental situations, the light field is our only source of informa- 
tion, while the system dynamics must be deduced indirectly from the respective measurement 
protocols. Furthermore, the type of measurement we choose to perform on the photons may 
influence the dynamics of the system. Therefore, a closer look at the processes which allow 



us to measure the emitted photons may be helpful. For this purpose we should start with the 
complete Schroedinger equation of the system and the field. The evolution of this state will 
then be analyzed with respect to a representation adopted to given measurement scenarios. 

II. BASIC MODEL 

In the most simple case, we have a localized two level system, consisting of the excited 
state \E) and the ground state \G) with a transition frequency ouq, and a one dimensional 
field with linear dispersion u = c\k\ coupled by a local interaction. 

The field can be separated into two distinct branches, one with group velocity c for 
k > and one with group velocity -c for k < 0. In the following, we will further simplify the 
problem by considering only the branch with positive group velocity, extending it to negative 
values of k and uj as shown in figure (1). The introduction of negative frequencies can be 
justified on the ground that, since the emission of photons with zero or negative energy 
would violate energy conservation, such processes will only contribute to the dynamics on 
time scales smaller than 1/ujq. Therefore, the amplitudes of those spurious field states should 
be small compared to those near uj = uj ||. 

The relevant states for the emission process are the product state of the excited system 
state and the field vacuum, \E;vac), and the product states of the ground state and the 1 
photon field states \G; k). Using these states as basis, the Hamiltonian of our model is 

H = hujQ | E; vac. ) (E; vac. \ 

/+oo 
dkhkc\G; k)(G;k\ 
-oo 

/+oo 
dkh[g\G;k)(E;vac.\ + g*\E;vac.)(G; k\] (1) 
-oo 

This Hamiltonian is a simple version of the Hamiltonian used in Wigner-Weisskopf theory. 
Therefore, we can proceed as in |4]]. The time dependent Schroedinger equation is 

— (E;vac.\vjj(t)) = —icu (E;vac.\ip(t)) 
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/+oo 
dk(G;k\m) (2) 
-oo 

±(G;k\m) = -i\k\c(G;k\m) 

- ig{E-vac.\^{t)) (3) 

If we choose the initial condition \ip(t = 0)) = \E;vac), the result of the integrated 
Schroedinger equation for the one photon part is 

(G; A#(t)> = -ig t dt' e'^'^ (E; vac. U(t')) (4) 
Jo 

This result can be substituted into the equation for the amplitude of the excited state. By 
using b(t) := e t<JJ ° t (E;vac.\'4)(t)) one obtains 



d_ 

dt 



POO /"+00 ( w M 

b(t) = -\g\ 2 dt'b(t') dke-^^W-V ( 5 ) 

JO J-oo 



The integral over k, again, results in a delta function, and the evolution becomes local in 
time: 



(E-vac.\ij(t)) =e ( - l ^' r/2)t (7) 
The amplitudes calculated for the 1 photon states of our model are: 

1 _ p -Tt/2 i(kc-w )t 

We can test our assumption that the amplitudes for k < are much smaller than those 
near k = u /c by calculating the k-space probability densities: 

1/^,1,^x12 i .2 1 - 2cos(kct - u t)e~ rt/2 + e" r * 

KG;t| * W>l =l91 + (9) 

For t >> 1/r, this is the familiar Lorentzian distribution. Thus, for ojq » T, the contri- 
bution of field states with u < will be negligible. 
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On smaller timescales, there are contributions from k = — oo to k — +00 due to the 
energy-time uncertainty. In a model with a more realistic dispersion relation, the k=0 
component have a groupvelocity of 0. Consequently, there will be reabsorption of such com- 
ponents in the real space formulation, which leads to an interaction of the system with itself. 
The main effect of this interaction would be an energy shift in the spectrum, comparabel 
to the lamb shift in quantumelectrodynamics. However, since our model does not include 
such processes, there is no such energy shift and the Lorentzian distribution of the emission 
centers on c^o- 

It is now possible to calculate the real space amplitudes by Fouriertransforming the k 
space result: 

-i x fL e {T/2+iu Q ){x/c-t) for Q <x <ct 



(G;x\m) = { 



(10) 

otherwise 



Figure (2) shows the spatiotemporal distribution of the field state. This function can be 
interpreted as the real space probability density for photon detection. 



III. APPLICATIONS 

A. Extension to 3 Dimensions 

In order to apply this result to realistic conditions in photon spectroscopy, it is necessary 
to give a more detailed interpretation of the \k) states. In typical cases, the emission will be 
into an unrestricted three dimensional light field. Therefore, an infinite number of degenerate 
modes with the same absolute value of k is available. In the plane wave basis with wavevector 
k = |k|k and a linear polarization given by e p , the interaction part of the Hamiltonian is 

Him = I d\k\dkJ2 ep hg(k,e p )\G;k,e p )(E;vac.\ 

+ hg*{k,e p )\E;vac.)(G;k,e p \ (11) 

However, instead of using this plane wave basis, we may transform to another basis in which 
only a single mode per energy level interacts with the system. This new mode can be 
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determined directly from the interaction part of the Hamiltonian: 



int 



/ G?|k|7k?|G; k)(E;vac.\ + hg\E; vac.)(G; k\ 



(12) 



with \G;k) = /rfkEe p ^|G;k,e p > 
and # 2 = /dk£ Gp |#(k,e p )| 2 



(13) 



(14) 



In many situations, it may suffice to inspect the symmetry of the system and of the field. 
For example, the spherical symmetry of an atomic system allows us to classify both the 
system and the field states according to the quantum numbers of the angular momentum, 1 
and m (for details, see J7|, ||). In the following, we restrict the atomic states to \S) (1=0) 
and \P) (1=1). A dipole transition emits only photons of Z = 1 with m = 0, ±1. Since the 
total angular momentum is conserved, the system must change its quantum numbers 1 and 
m accordingly. 

Since for radii much larger than the wavelength spherical waves approach plane waves, 
the previous interpretation of the wave function in real space still applies, except in the 
immediate vicinity of the object. Replacing x by r, the |r) states now describe a photon 
at a distance r from the system, while the angular dependence and the polarization of the 
photon is given by the photon state quantum number m. With the z-axis as quantization 
axis, the amplitudes are 



The condition for transversality of the light field, e p k = 0, yields the angular dependence of 
the emitted intensity characteristic for dipole radiation. 



A simple application of this model is the standard quantum beat scenario (three level 
system) in which an atom in a magnetic field is excited from the groundstate \G) = \S) 




(15) 



(16) 



B. Three Level System 
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with a short laserpulse polarized in a direction orthogonal to the magnetic field. The atom 
is left in a superposition of the m — +1 and the m — — 1 sub levels of the excited P state 
\E±) = \P, m — ±1) . Since we can apply our theory to both the Am = +1 and the Am = — 1 
transitions separately and since we are dealing with the linear dynamics of the Schroedinger 
equation, we can immediately write down the complete evolution of both the field and the 
system in local representation: 

\iP(t = 0)) = -={\P,7n= +1; vac.) + \P,m = -1; vac.)) (17) 
v2 

\iP(t)) = -^ r e- r+t/2 ~ luJ+t \P,m = +l;vac.) 
v2 

+ -^c~ T - t/2 - iuj - t \P,m = — 1; vac.) 

ret Fp 



JL f C ( /i± e (r + /2+ia, + )(r/ C -*)| (S - = +1 \ 

v2 jo V c 



V2 

+ ^eF-/ 2+i "-X r / c -V\S-,r,m= -l))dr (18) 

Written in the m basis of the light field and the m basis of the system, quantum beats do 
not appear in the evolution of either subsystem. A simple transformation to another basis 
in atom-, 

\P X ) = ^=(\P,m= +1) + \P,m= -1)) (19) 
\P y ) = A^{\P,m= +1) - \P,m= -1)) (20) 

and field-states 

\r,d x ) = —f=(\r,m = +1) + \r,m= -1)) (21) 

\r, d y ) = -J=(|r, m = +1) - \r, m = -1)) (22) 

shows, however, that beats do occur then, contrary to the statement made in 0, that 
orthogonal dipoles should definitely exclude beats. The wave function written in this basis 
for the case of T + = T_ := T now reads: 
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\^(t)) = e- rt/2 ' iidt (cos(5Lot)\P x ;vac.} + sin(5ut)\P y ;vac.)) 
~ *{^r e(r/2+ " )(r/C_ * ) (cos(^(r/c - t))\S; r, d x ) 
+ sin(5uj{r/c-t))\S;r,d y ))dr (23) 

where u> = and <5c<j = , This representation of the evolution clearly shows beats 

in both the system and the field. Figure (3) shows the real space probability distribution 
of this linearly polarized photon. Of course, whether beats are observed or not depends on 
the type of measurement performed. 



C. System-Field Correlations 

Another application of this model is the description of entangled states between field 

and system if the decay is from a single excited level \E) = \S) into two (or possibly 

more) alternative groundstates \G±) = \P,m — ±1). Although the situation may at first 

appear very similar to the quantum beat scenario, we cannot single out two paths, since the 

initial state cannot be separated into two components, each belonging to only one path. A 

possibility is to view the decay as a single path, leading from the excited state to a linear 

combination of the two groundstates, —7== IP, m = +1) H — , r ~ „ IP, m = — 1). However, 

V r ++ r - 1 ' ' ^r 2 + +r 2 _ 1 ' ' ' 

the two groundstates will generally not be degenerate. Therefore, the state of the system 
will start to evolve in time at the instant of the emission, t-r/c, causing a dependence of the 
system state on the field state |r). The measurement of a photon at r implies that a time of 
r/c has elapsed since emission, so that the state of the total system would be 

\G; r) := 1 (T + e luJ+r/c \P, m = +l;r,m = -1) 



+ r_e iW - r/c |P, m = -l;r,m = +l)) (24) 

Thus we can use the fact that our model treats emissions as instantaneous to replace the 
ground state with the temporal evolution of this two level subsystem. 
With T = T + + T_ we obtain the result of the complete evolution as 
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= e - rt/2 - iwot \E;vac.) -i 
= e -F+/*+r-/M*>)t\S ]V ac) 



r r ct 



c Jo 



jr/2+iw )(r/c-t)\Q.^ r ^ ^ r 



T,+r. 



ct 



3 (r + /2+r_/2+^ )(r/c-t)_ 



ri + r? 



{T + e iuJ+{r/c - t) \P,m = +l;r,m 



-1) 



+ T^e i "- {r/c - t) \P,m = -l;r,m = +1)) dr (25) 

The one photon part of this wave function carries a strong system-field correlation. In the 
case r + = T_ = T/2, any linear combination of \P, m = +1) and \P, m = —1) states 
may be observed, depending on the measurement results in the field. For example, we may 
transform to the field basis \r, d x ), \r, d y ) and the atomic basis \P X ), \P y ) as in section 3.2.: 

\^(t)} = e- {r/2+luJo)t \S;vac.) 

-ij- [** e F/2+uj +iu>)(r/c-t) J_ 



C Jo y/2 

(cos(6(jj(r/c-t))(\P x ;r, d x ) + \P y ;r, d y )) 

+ sm(6uj(r/c - t))(\P x ; r, d y ) - \P y ; r, d x ») dr (26) 

Again, we use the notation a) = aJ+ ^ a; ~ and 5u; = u '+~ u} - . While the total probability 
of measuring a photon with d x polarization does not oscillate, those photon detections 
correlated with an atomic state of \P X ) and \P y ) both display beats as shown in Figure (4). 
The beats could therefore be seen by measuring the atomic state and separating out the 
photon detections accordingly. This seems to be at variance with the expectation that in 
this type of three level system, beats should be absent because a measurement of the atomic 
state would reveal the decay channel, u + or w_ 0. However, measuring the \P X ) and \P y ) 
states of the atomic system reveals only the phase of the beats, preventing the determination 
of the decay channel. Coincidence measurements of this type can therefore be considered 
as quantum eraser measurements as described in fH| and [0 for the case of 2 spatially 
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separated sources of photonscattering , since the measurement of \P X ) and \P y ) effectively 
erases the information to be gained by measuring \P, m = +1) and \P,m = —1). 

The two cases described above and in section 3.2, respectively, may be combined into a 
single four level system with a cascade decay leading from an excited level through two (or 
more) intermediate levels to the ground state. This model could thus be used to describe 
the quantum beats observed in cascade decays fll2| . The resulting wave function will include 
a two photon part, so it may be necessary to consider the boson nature of the photons when 
transforming to another measurement base. 

Since the emission process is local in time and there is no interaction between the pho- 
tons, our present model can be adapted to any other scenario by simply adding the effects 
of the possible channels of decay. In this way, it should be possible to develop a better 
understanding of the role which the optical measurement apparatus plays in quantum mea- 
surements. 



IV. MEASUREMENT PROCESSES 

For actual measurements the atomic object and the field is embedded into a dissipative 
environment, usually a filter and a detector, which, again, may be specified by its interaction 
with the light field. The type of measurement process can be included in the form of a 
projection acting on the field part of the wavefunction. The choice of base states used 
should thus be made with a special experimental setup in mind. For example, a polarizer- 
detector setup can be considered as a projective measurement on |r; d x ) states, where r is 
the distance between the system and the detector. With the 3-level system described by 
equation (23), this would allow us to observe quantum beats. A more realistic representation 
would consider the limited area covered by the detector, meaning that actually a state of 
the form (^ s /a : \r,d x ) + \/l — a\r, I > 1)) is measured, where a is the proportion of photons 
emitted into the angle covered by the detector. \r,l > 1) represents a superposition of 
multipole states with I > 1. It ensures, that the measured state exists only in a limited area 
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by interfering with \r, d x ) in such a way, that real space components outside the emission 
angles covered by the detector vanish. 

In our model, the field dynamics in real space rigidly move the field states away from 
the system at the speed of light (see Fig.l). Therefore, the state of the field represents a 
temporal record of the field dynamics at the system. It does not matter at what distance r 
we choose to measure, since a shift in r is equivalent to a delay time At = r/c. If we are only 
interested in time resolved spectroscopy, we may simulate the photon measurement directly 
at the system. This corresponds to the measurement approaches of and 0. 

However, as we have calculated the whole wavefunction of the field, we can also con- 
sider frequency resolved measurements. Most frequency selective filters are based on an 
interference between optical paths of different length. This corresponds to a linear com- 
bination of different |r) states of the field. A Michelson interferometer would produce a 
measurement of ^(l r i) + l r 2)) in the detector, while a Fabry- Perot interferometer measures 
^ 1 l i? 4 Y^=i R 2n \r = nd), R being the reflectivity and d the distance between the panels. 

In the case of cascades, the consecutive decays will produce a many photon wavefunc- 
tion. Nevertheless, we can still apply our theory by using the many photon real space base 
|ri,r 2 , ...). Coincidence measurements such as applied in |T2| can now be represented by a 
measurement base of |ro,ro + cr), with r as the delay time. 



V. CONCLUSIONS 

We have presented an exactly solvable model for the unitary temporal evolution of the 
total system-field state. The assumptions used in the model are equivalent to the approxima- 
tions of Wigner-Weisskopf theory. However, our approach enables us to provide a physical 
interpretation of these approximations in terms of locality and constant group velocity. 

Also, it is possible to visualize the evolution of the wave function in space and time, 
allowing us to take a closer look both at time resolved spectroscopy and at the effect of 
interferometric filters. This clearly reveals the special role of the measurement process in 
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quantum mechanics, which is often concealed by the application of approximations such as 
the bath approximation or the omission of correlations. 

The possibility of adapting the model to different situations may provide an opportunity 
to study the emission process in complex quantum systems such as molecules and nanos- 
tructures. 
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FIGURES 



Fig. 1. Extension of the linear dispersion uj = \k\c (solid line) to negative Frequencies (dashed 
line). The dotted line marks the resonant frequency of the system. 

Fig. 2. Spatial distribution of the probability of photon detection for t = In2/Y (solid line), 
t = 2ln2/T (dashed line) and t = 3ln2/T (dotted line). 

Fig. 3. Spatial distribution of the probability of photon detection for lineary polarized photons 
in a quantum beat experiment. 

Fig. 4. Beats correlated with the \P X ) state (solid line) and the \P y ) state (dashed line) of 
the atomic system. The dotted line shows the sum of both probabilitydensities, which displays no 
beats. 
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r in units of In2*c/T 



